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Most cosmic-ray air shower arrays have focused on detecting electromagnetic shower particles and low energy 
muons. A few groups (most notably MACRO + EASTOP and SPASE + AMANDA) have studied the high 
energy muon component of showers. However, these experiments had small solid angles, and did not study muons 
far from the core. The IceTop + IceCube combination, with its 1 km 2 muon detection area can study muons far 
from the shower core. IceCube can measure their energy loss (dE/dx), and hence their energy. With the energy, 
and the known distribution of production heights, the transverse momentum (pr) spectrum of high pr muons can 
be determined. The production of these muons is calculable in perturbative QCD, so the measured muon spectra 
can be used to probe the composition of incident cosmic-rays. 



1. Introduction 

Despite much effort, the composition of high- 
energy (above 10 15 eV) cosmic-rays is poorly 
known. Satellite and balloon experiments have 
insufficient exposure for high-statistics measure- 
ments. Terrestrial experiments have used the 
electromagnetic and low energy muon compo- 
nents of air showers for composition measure- 
ments. The ratio of electromagnetic to muon en- 
ergy may be sensitive to the cosmic-ray composi- 
tion, but it is also sensitive to the physics models 
used to simulate the air shower. One critical com- 
ponent of these models is the forward production 
of hadrons in high-energy interactions. Most of 
these hadrons are produced at low transverse mo- 
mentum (pt), where perturbative QCD (pQCD) 
is not applicable. A variety of non-pertubative 
calculations and extrapolations are used to model 
low px particle production, and thereby make 
predictions about cosmic-ray composition. 

Here, we present a complementary approach 
to study air showers, using high energy, high px 
muons to study cosmic-ray composition. The cal- 
culations use pQCD, and so should be well un- 
derstood. High px muons are produced predom- 
inantly in the initial cosmic-ray interaction pQ, 
from semileptonic decays of heavy quarks, and 
from decays of high pt pions and kaons produced 
in jet fragmentation. The rates for muons with 
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Pt above a few GeV is calculable in pQCD; the 
spectrum depends on the parton composition of 
the incident ions. For a fixed shower energy, the 
energy per nucleon drops as the atomic number 
rises, substantially altering the parton density, 
and so changing the muon spectrum. This com- 
position may be inferred from the high pt muon 
spectrum. The low— x parton distributions in ni- 
trogen also contribute in some kinematic areas; 
this may be an additional study topic. 

2. Experimental Technique 

The study of high px muons requires a surface 
air shower array combined with a large under- 
ground muon detector. The surface array mea- 
sures the shower energy, core position and in- 
cident direction, and the underground detector 
measures the energy and position of high-energy 
muons. Previous experiments have studied high 
energy muons in air showers, but, with relatively 
small underground detectors [2] [3] ■ These exper- 
iments did not make use of the distance between 
the muon and the air shower core. 

IceCube and IceTop comprise a 1 km 2 sur- 
face air shower array and a 1 km 3 muon detector 
[4]. The muon detector is big enough to observe 
muons far from the shower core. Together, the 
combination can determine the key elements of 
the event. 

IceTop will measure the shower energy, core po- 



2 



S. Klein 



sition, and arrival direction. It will consist of 160 
ice-filled tanks spread over a 1 km 2 area[5j. It has 
an energy threshold of about 300 TeV. 

IceCube will consist of up to 4800 optical mod- 
ules in 1 km 3 of ice, at depths from 1450 to 2450 
meters. The combined acceptance will be about 
0.3 km 2 sr [5]. IceCube will measure the en- 
ergy, position and direction of muons. For ver- 
tical muons, the energy threshold is about 500 
GeV Muon energy, is measured by determin- 
ing the muon specific energy loss (dE/dx). For 
Ef, > 1 TeV, dE/dx scales with E^. 

In a sufficiently high energy event, IceCube will 
observe multiple muons. Most of these muons 
will have low px, and so will cluster around the 
shower core. In this high-density region, it may 
not be possible to reconstruct individual muons; 
the bundle will be reconstructed as a single light 
source. However, far from the core, where the 
muon density is lower, it should be possible to 
reconstruct individual tracks. The separation re- 
quired to resolve individual tracks is unknown, 
but 100 meters (comparable to the spacing be- 
tween optical module strings) seems like a safe 
value. It is more than 3 times the effective light 
scattering length of about 30 m [BJ. 

With the muon energy and distance from the 
core determined, the muon pt may be calculated: 

Pt = — — (1) 

where E^ is the muon energy, d c is its distance 
from the core, and h is the distance from IceCube 
to the site of the initial cosmic-ray interaction in 
the atmosphere. The value of h is not determined 
on an event-by-event basis, but its average value 
(w 30 km), zenith angle dependence and distribu- 
tion are well known; the event-to-event variation 
and slight composition dependence can be con- 
sidered in determining the pr spectrum. 

With the conservative d c > 100 m requirement 
and a 1 TeV muon, IceCube can study the spec- 
trum for pt > 3 GeV/c, covering most of the 
interesting high pt region. As increases, so 
does the minimum accessible px] for E„ = 10 
TeV, p T > 30 GeV/c, while for a 50 TeV muon, 
Pt > 150 GeV/c. Since higher energy muons 
produce more Cherenkov light, they are easier 



Table 1 

Muon rates for the TIG95 calculation, for 1 year 
(3 x 10 7 s) with 0.3 km 2 sr acceptance. 



Energy 


Prompt 


Non-Prompt 


Threshold 


Rate 


Rate 


10 15 eV 


1.5 


0.5 


10 14 eV 


1,000 


18,500 


10 13 eV 


56,000 


10 7 


10 12 eV 


600,000 


7 x 10 8 



to track, and it is likely that, for higher energy 
muons, a smaller d c cut could be used. 

The systematic errors in these measurements 
remain to be studied. Some components are: 

1) Uncertainty in the absolute cosmic-ray flux 
and energy scale. The flux factors out, since we 
will only use events observed by IceTop. The en- 
ergy scale introduces an uncertainty in the scale 
of x measurements. 

2) Error on the core position and extrapolation 
to depth. For muons far from the core, the frac- 
tional error is small. In the first IceCube string, 
the systematic offset between IceCube and IceTop 
is less than 1 degree [7]. 

3) Uncertainty in the muon energy and posi- 
tion, due to stochastic interactions, multiple scat- 
tering and other factors. Multiple scattering con- 
tributes to d c , but, far from the core, d c is domi- 
nated by the muon pt- 

3. Rates 

Thunman, Ingelman and Gondolo (TIG) calcu- 
lated the prompt (charm only) and non-prompt 
muon production using a PYTHIA Monte Carlo 
simulation, with the MRS-G parton distribution 
functions, leading-order pQCD cross-sections for 
qq — » QQ and gg — » QQ, and standard charm 
quark hadronization and decay models [8]. Par- 
ticles were propagated in the atmosphere using 
transport equations. Table Q] shows the expected 
muon rate in the combined IceCube -IceTop ac- 
ceptance for different energy thresholds. Calcu- 
lations, with a fixed E~ 3 7 energy spectrum find 
more non-prompt muons at low energies [S]. 

A newer calculation uses next-to-leading-order 
pQCD calculations and the CTEQ3 parton dis- 
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tributions. It finds prompt rates that are compa- 
rable to the earlier calculation at 1 TeV, but are 
higher at higher energies; at 100 TeV, the prompt 
cross section is about 8 times higher [T] . The dif- 
ference is due largely to the different low— a; be- 
havior of the two parton distributions. 

Bottom quark production was not included. 
Calculations for the LHC (at a comparable en- 
ergy) find that, for muon px > 2 GeV/c, bottom 
contributes a larger muon signal than charm 1 . 

For 2£ M > 50 TeV, the accompanying shower 
should almost always trigger IceTop; at lower 
muon energies, some of the showers may be too 
small to trigger IceTop, so will not be seen. 

The prompt signal is significant for < 10 14 
eV, although smaller than the non-prompt signal. 
A px cut should eliminate all of the soft (non- 
perturbative) non-prompt signal, leaving muons 
from heavy quarks and high px ft and K. Based 
on calculations for the LHC, a cut on pr > 2 — 4 
GeV/c will also eliminate 90-98% of the prompt 
muons [TO], Although a drastic reduction, this 
still leaves an interesting sample. The final cut 
on px (or d c ) will depend on the detector 2-track 
(core + distant muon) separation capability. 

4. Muon spectrum Analysis 

Measurement of the muon pt spectra has some 
similarities with the lepton spectra studies done 
at the Relativistic Heavy Ion Collider (RHIC). 
The pt spectrum of leptons produced in proton- 
proton, deuteron-gold and heavy-ion collisions 
have been studied. 

The RHIC experiments fit their data using a 
multi-component 'cocktail.' For electrons, the 
cocktail consists of ir° and 77 Dalitz decays, 
7 — > e + e~ conversions, leptonic decays of vector 
mesons, plus semileptonic decays of heavy mesons 
and baryons [11] . For muons, the cocktail consists 
of 7r, K and heavy quark decays |12j . For muons, 
the 7r and K decay fraction is reduced with vertex 
cuts, sample, so this result is not directly relevant 
to IceCube. However, the electron analysis seems 
quite relevant. The fraction of prompt electrons 
rises with the electron pt; for pt > 5 GeV/c, 
prompt electrons are dominant. Additional confi- 
dence in these studies comes from the good agree- 
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Figure 1. Parton distribution for a 10 17 eV cos- 
mic ray for quarks (top), gluons (middle) and an- 
tiquarks (bottom). The solid lines are for hydro- 
gen (A — 1), while the dashed lines are for A = 10 
nuclei. These curves are based on the MRST99 
parton distributions [14] at Q 2 = 1000 GeV 2 . 



ment seen between high pt ft data and pQCD 
calculations [13]. 

Simple arguments predict this dominance. 
Light meson production is predominantly soft; 
da/dpT falls exponentially with (pt) ~ m, m be- 
ing the meson mass. In contrast, heavy quark 
production is described by pQCD, which gives a 
power law pt spectrum; at high enough pt, this 
will dominate over any exponential. pQCD pro- 
cesses also produce high pt ft and K, but these 
mesons are only a small fraction of the total pro- 
duction. 

A similar approach could apply to the pt spec- 
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trum of muons from air showers, with a cocktail 
of perturbative and non-perturbative non-prompt 
muons, and prompt muons. For px above a few 
GeV/c, the non-perturbative component will be 
gone, and the remaining events could be fit to a 
power law spectrum. Spectra could be measured 
for different muon energies and shower energies. 
The muon energies are related to the muon ra- 
pidity in the center-of-mass frame. 



5. Composition Determination 

Fig. 1 compares the parton energy spectrum 
(in the target frame) for a 10 17 eV cosmic ray, 
for hydrogen (A = 1) and an A = 10 nucleus. 
The maximum parton energy scales as 1/A; this 
determines the maximum parton-parton center of 
mass energy. 

For a given collision, the x values of the parton 
densities that are probed depend on the kinemat- 
ics of the produced partons. So, different muon 
energies and px are sensitive to different x re- 
gions. 

In any collision, the maximum possible muon 
energy is the maximum incident parton energy. 
The maximum pj< is half the parton-parton 
center-of mass energy, W = yj2Expx^m v . Here 
E is cosmic-ray energy, xt and xj\r are the x of 
projectile and target partons, and m p is the mass 
of the proton. The corresponding spectra are de- 
termined by a calculation that includes the kine- 
matics of the parton production, fragmentation, 
and semileptonic decays. 

It may be worth giving a few examples of the 
range of x values probed, using a a grossly sim- 
plified model where the muon takes half of the 
energy and pt of the parton produced in the col- 
lisions. For an incident 10 18 eV proton produc- 
ing a 10 16 eV muon with p T = 1 GeV (this muon 
would be in the core, and only distinguishable due 
to its huge dE/dx), xp = 0.01, x^ = 1.5 x 10 -6 
and Q 2 « 30 GeV 2 . With the same incident par- 
ticle and muon energy, a muon with pt = 100 
GeV would come from a collision with xp = 10~ 3 , 
x N = 8 x 10~ 2 and Q 2 = 10 5 GeV 2 . 



6. Conclusions 

Studies of high px muon production in cosmic- 
ray air showers appears feasible with IceCube and 
IceTop combined. IceTop can measure the air 
shower energy, incident direction, and core po- 
sition, while IceCube will measure muon energy 
and distance from the shower core; these data can 
be used to determine the muon px- For px above 
a few GeV, the muon px spectrum can be inter- 
preted in terms of perturbative QCD plus and 
fragmentation functions. The muon px spectrum 
should be sensitive to the composition of incident 
cosmic rays. 

It is a pleasure to acknowledge useful comments 
from Xinhua Bai and Tom Gaisser and Teresa 
Montaruli. 
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